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Abstract 

We study Brownian particle motion in a double-well potential driven by an ac force. This system exhibits the 
phenomenon of stochastic resonance. Distribution of work done on the system over a drive period in the time asymptotic 
regime have been calculated. We show that fluctuations in the input energy or work done dominate the mean value. 
The mean value of work done over a period as a function of noise strength can also be used to characterise stochastic 
resonance in the system. We also discuss the validity of steady state fluctuation theorems in this particular system. 
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1. Introduction : 

Stochastic resonance (SR) refers to the enhanced 
response of a system to a small deterministic peri- 
odic forcing in the presence of an optimal amount 
of noise [Ij . It often occurs in bistable and excitable 
systems with subthreshold inputs. Noise plays a con- 
structive role in this phenomenon. This is due to a co- 
operative interplay between nonlincarity of the sys- 
tem, input signal and noise. Thus the power from the 
whole spectrum is pumped into a single mode that is 
coherent with the external driving force. Because of 
its generic nature, this phenomenon, boasts applica- 
tions in almost all areas of natural sciences. Differ- 
ent quantifiers of SR have been discussed [3i|4l|5|6f7] 
as regards to its validity as a bonafide resonance, 
i.e., the resonance phenomena as a function of noise 
strength as well as the driving frequency. Hysteresis 
loop area[3], input energy or the work done on the 
system per cvcle [4l5| and area under the first peak in 
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the residence time distribution [2| have turned out to 
be good quantifiers characterizing SR as a bonafide 
resonance. Recently it has been shown that the dif- 
ferent quantifiers of SR are mathematically related 
to each other fB^. Earlier studies on the work done 
in a periodically driven bistable system by an ex- 
ternal agent has established that the average work 
peaks around the resonance, as expected [4|5| . In 
our present work, we study the fluctuations and na- 
ture of the probability distributions of work done 
over a period in a driven double-well system in the 
time asymptotic regime (where all averaged c^uan- 
tities become periodic in time with a period of the 
external drive). 

Motivation for this study also comes from re- 
cent interest in the so called fluctuation theorems 
and transient violation of second law in systems 
driven to a nonequilibrium state |8I9"] . Fluctuation 
theorems describe properties of the distribution of 
various nonequilibrium quantities, such as work, 
heat, entropy, etc. [8 9 10 . It is important to note 
that these theorems are rigorous and applicable 
to systems driven arbitrarily far away from equi- 
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librium. These theorems can be extremely useful 
in analysing the role of fluctuations on the perfor- 
mance characteristics of engines at nanoscales (e.g., 
molecular motors). In these tiny systems interac- 
tions with the environment are dominated by ther- 
mal fluctuations |8j. Moreover fluctuations in phys- 
ical quantities are more than the mean values and 
large fluctuations may occur which can lead to un- 
expected consequences. It has been shown recently 
that overdamped driven harmonic oscillator (linear 
system) satisfies a steady state fluctuation theorem 
(SSFT) for the work done over a single period {Wp) 
in the time asymptotic regime[TT|. SSFT implies 
that the probability distribution P{Wp) of work 
Wp satisfies the following relation [10|ll|12|13|14j . 
namely, 



PjWp) 
P{-Wp) 



(1) 



We show that this is an artifact of the linearity of 
the model and in general does not hold for driven 
nonlinear systems. However, we show that SSFT for 
the work is indeed satisfied, if one instead considers 
the work done over a large number of cycles [M]. 
Our study of input energy fluctuations also reveals 
that the relative variance is larger than one across 
SR indicating that fluctuations dominate the mean 
value. 

2. The Model: 



The overdamped Langevin equation of a particle 
in a double-well potential in the presence of a time 
periodic force is given bv|12j 



dx 

It ~ ' 
where 



dUjx) 
dx 



(2) 



U{x) 



Ax sin Lut. 



(3) 



X X 

T " y 

We have set friction coefficient (7) to unity. The 
random force field ^ (t) is a zero mean Gaussian white 
noise, i.e, 



mat )) = DoS{t-t). 

where Dq = 2'jkBT is the 
static double- well potential V{x) 



(4) 

where Dq = 2'jkBT is the strength of the noise. The 
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Fig. 1. Figure shows the variation of average work Wp done 
per period with temperature D. Inset shows the correspond- 
ing plot of relative variance Vr- The parameters are A=0.1, 
and LU = 0.1. 

a barrier height AV — 0.25 between the two sym- 
metrical wells (minima) located at distances Xm = 
±1. We consider the case of weak forcing < 
Ay. All the physical quantities are taken in dimen- 
sionless units as prescribed in references |4I5| . 



The work done by the external drive on the system 
or the input energy fr^ = — j is defined as [4l5j 



" I dt 

to 



^—Aijj j x{t) COS ujtdt 
to 



(5) 



This follows from the stochastic energetics formal- 
ism developed by Sekimotojl6j. 

Numerical simulation of this model was carried 
out by using Huen's method [T7j. To calculate the 
work done we first evolve the system and neglect 
initial transients and then work done over a period 
[Wp) is calculated (eqn.(5)). The values of Wp are 
different for different cycles. To get better statistics 
we have calculated Wp for 100000 different (some- 
times even more) cycles. 
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Fig. 2. Probability distribution for various values of D, other parameters being the same as in Fig.l. 



3. Results and Discussion: 

In Fig.l we plot the average work done (Wp) over 
a single period (in the time asymptotic regime) as a 
fmiction of noise strength D(= kT = i, in dimen- 
sionless units) for low amplitude driving of strength 
A ~ 0.1 and for frequency oj — 0.1. We clearly see 
that input energy exhibits a peak as a function of 
D. This can be attributed to the synchronized es- 
cape from the potential well with the external pe- 
riodic drive. We reproduce exactly the same figure 
as in reference[5], which has been obtained using a 
different numerical method. 

For the case where noise is small the particle 
remains in one of the potential wells for a longer 
time(intrawell motion dominates) with an occas- 
sional random jump to another well (Kramers' es- 
cape over the barrier) as a function of time. When 
the noise is very strong a large number of interwell 
random switches occur for each period of the si- 
nusoid and the systems response is again random 
(interwell motion dominates). In between these 
two conditions, surprisingly, there exists an opti- 
mal value of the noise that cooperatively concurs 
with the periodic forcing to make almost exactly 
one switch per period. Hence interwell motion will 
be in synchrony with the input signal. Quantita- 
tively this condition is determined by the matching 
of the two time scales [115] . namely, the period of 
the input signal and the Kramers' escape rate r^, 



1 



OJ 

- = rk: 



1 

with ri- — — 
2n 



At 



the optimal value of noise the input energy is max- 



imum. The nature of this maximum is a function 
of other physical parameters as discussed in detail 
in|4|5]. In the inset of Fig.l we have plotted the rel- 



ative variance Vr = 



{w/)-{Wpy 



of work done 



Wp over a period as a function of noise strength. 

We see that Vr > 1 in the parameter range we 
have chosen. This shows that fluctuations dominate 
the mean value. Vr exhibits a minimum near the 
value of D at which SR occurs. It may be noted that 
when Vr becomes much larger than one the average 
value ceases to be a good physical variable (non self- 
averaging quantity) and hence to characterise the 
behavior of work one has to resort to the study of 
the entire distribution function. In the time periodic 
state the average work done is dissipated into the 
system as heat. Thus one can identify (Wp) as a hys- 
teresis loss in the medium. However, it may be noted 
that fluctuations of work done cannot be identified 
with heat fluctuations [llll2j . 

Now we will turn to the understanding of the na- 
ture of probabihty distribution P{Wp) ofWp for var- 
ious values of noise strength D across SR peak. For 
this we have plotted P{Wp) as a function of Wpin 
Figs. 2a and 2b for various values of D. For low values 
of D (^ 0.02) particle dynamics is mostly confined 
to a small amplitude intrawell oscillatory motion. 
Hence the distribution is closer to the Gaussian. For 
slightly higer values of D(~ 0.04 and 0.06) the par- 
ticle makes occassional random excursions into the 
other well as a function of time. This is clearly re- 
flected as a small asymmetry and a hump (interwell 
motion) in the plot of P{Wp). It is interesting to 
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Fig. 3. 

Plot of P{Wp) (line with circles), 
P{~Wp)exp{(3Wp){\ine with crosses) vs. Wp. Plots 
are for D=0.02(a), D=0.16(b) and D=0.5(c) re- 
spectively. For all plots w = 0.1 and A — 0.1. 




pWnp 



Fig. 4. Plot of ln{PiW„p)/P{-Wnp)} vs. l3W„p (over 10 
periods). Inset shows the corresponding plot of Wnp (over 5 
periods). The plots are for D=0.16,tj = 0.1 and A=0.1 

note that there is substantial weight towards nega- 
tive values of Wp . The negative values of Wp corre- 
sponds to transient second law violating hysteresis 
loops. As we increase D further the interwell dynam- 
ics starts playing a prominent role. It is interesting 
to note that position of the first and the second peak 
does not shift much, however, its width increases. 
The weight of the probability distribution towards 
higher negative values of Wp increases. On further 
increasing D (Fig. 2b) finally a single peak structure 
appears. The shape of the single peak structure is 
closer to a Gaussian at high D('-~^ 0.5). For such high 
values of D particle makes several random excur- 
sions between the two wells during a single time pe- 
riod of an external drive. Hence total work can be 
treated as an addition of independent increments of 
work. Then the central limit theorem leads us to ex- 
pect that the distribution of work will be approxi- 
mately Gaussian. This is indeed the case. The pa- 
rameter range for which we have explored the distri- 
butions are broad and the variance of fluctuations 
dominate over the average values (Fig.l, inset). This 
is one of our main observations. 

Finally we discuss the applicability of SSFT for 
the work done over a single period. For this we have 
plotted in Fig.3, P(Wp) and Pi-Wp)exp{pWp) for 
three different values of D. For small D (~ 0.02), 
Fig.3a and for large D (~ 0.5), Fig.3c PiWp) and 
P[—Wp)exp{(3Wp) are almost equal, thus satisfying 
SSFT for work distribution. For the intermediate 
value of D (-0.16), Fig.Sb we see that SSFT does 
not hold[18 . This clearly indicates that prediction 
of validity of SSFT for work done over a single period 
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is restricted only to a driven overdamped harmonic 
oscillator (linear problem) [11] . This will not be the 
case in the presence of nonlinearity. We consider a 
parameter regime where D — 0.6, lo — 0.01, A — Q.\ 
of Fig. 3, where SSFT over a single period does not 
hold. However, we calculate the work distribution 
Wnp for 5 and 10 periods. 

In Fig. 4 we have plotted ln{-p^^^^) as a func- 
tion of pWnp for periods 5 (inset) and 10. From 
the plots it is obvious that the slope is 1 and hence 
proving the validity of SSFT [14] . Moreover, the ob- 
served probability distribution is Gaussian and its 
variance V is related to the mean value (Wnp), i.e., 
V ^ ^{Wnp), which is a requisite for SSFT if the 
distribution is Gaussian (eqn. 24 offllj). 

Thus we see that work distribution satisfies SSFT 
provided we consider work over a large number of 
periods. The number of periods above which SSFT 
holds depends on the physical parameters. The finite 
time corrections are very complex|14] which we have 
not discussed here. 



4. Conclusion: 

We have studied a problem of work fluctuations in 
a driven double well system which exhibits SR. We 
have shown that across SR fluctuations of work cal- 
culated over a single period in the asymptotic time 
region dominates the mean value. The distributions 
exhibit significant weight towards negative values of 
work indicating transient second law violating tra- 
jectories. As opposed to the case of a linear problem 
of driven harmonic oscillators SSFT for work distri- 
bution is not satisfied for a work calculated over a 
single period. However, if the work is calculated over 
a large number of periods SSFT is indeed satisfied. 
We have also calculated the distributions of internal 
energy and dissipated heat over a cycle so as to ver- 
ify the extended form of SSFT [12]. The work will 
be reported elsewhere. 
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